function benchmark(D0i,lambda,gam,s_min,s_max,dist_s,phi,y1,y2,optfig)

% Benchmark with no runs/no DI

R1_grid = (1.001:0.001:2)';

J_FB    = zeros(size(R1_grid));
U_early = zeros(size(R1_grid));
U_late  = zeros(size(R1_grid));
for i = 1:length(R1_grid)
    R1 = R1_grid(i);
    [J_FB(i),U_early(i),U_late(i)] = fn_J_FB(R1,D0i,lambda,gam,s_min,s_max,dist_s,y1,y2,phi);
end

% Figure

% fig_first_best(R1_grid, J_FB, U_early, U_late, optfig)

[~,ind_max] = max(J_FB);

fprintf('Optimal benchmark first-best R1 is %6.3f \n\n', R1_grid(ind_max))

end
